Topological Inflation with Multiple Winding 



Andrew A. de Laix*, Mark Trodden^ and Tanmay Vachaspati* 

Particle Astrophysics Theory Group 
Department of Physics 
Case Western Reserve University 

10900 Euclid Avenue 
Cleveland, OH 44106-7079, USA. 



Abstract 



We analyze the core dynamics of critically coupled, superheavy gauge vor- 
tices in the (2+1) dimensional Einstein- Abelian Higgs system. By numerically 
solving the Einstein and field equations for various values of the symmetry 
breaking scale, we identify the regime in which static solutions cease to ex- 
ist and topological inflation begins. We explicitly include the topological 
winding of the vortices into the calculation and extract the dependence on 
the winding of the critical scale separating the static and inflating regimes. 
Extrapolation of our results suggests that topological inflation might occur 
within high winding strings formed at the Grand Unified scale. 
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I. INTRODUCTION 



Two of the most exciting possibilities suggested by modern particle cosmology are the 
production of topological defects and the existence of an inflationary stage in the early 
universe. In the last few years it has been pointed out independently by Vilenkin Jl] and 
by Linde || , that these two phenomena may be related. A necessary condition for inflation 
to occur is that the energy density of the universe be dominated by the vacuum energy 
of a homogeneous scalar field, called the "inflaton". Correspondingly, a central feature of 
topological defects is that their cores are regions of spacetime in which a scalar field is 
forced to sit out of its vacuum manifold, taking a value corresponding to a false vacuum. 
Hence, the cores of topological defects are regions in which energy density is trapped in the 
form of the vacuum energy of a scalar field. Vilenkin and Linde realized that if the energy 
density trapped in a defect core and the core radius itself were large enough, this energy 
could be considered uniform (ie. horizon sized), and the core would satisfy the conditions 
for inflation. This scenario is what is known as topological inflation. The advantage of this 
implementation is that, while traditionally the necessary conditions require the fine tuning 
of a scalar potential, here inflation is inevitable if the vacuum manifold of the theory at 
high temperatures satisfies a topological constraint. Thus, the question of initial conditions 
becomes one of topology. 

If defects are produced when a scalar field gets a vacuum expectation value equal to r\ 
in the early universe, then an approximate criterion proposed by Vilenkin and Linde for 
topological inflation to occur is rj > m p , where m p is the Planck mass. (Note that the 
energy density in the scalar field is proportional to At/ 4 which is assumed to be less than 
the Planck energy density.) The supposition that, under this condition, time-dependent 
defect solutions are inevitable, is supported by the traditional solutions for the spacetime 
around, for example, a cosmic string. At symmetry breaking scales significantly below the 
Planck scale, the spacetime around a static cosmic string is conical 0, with a deficit angle 
proportional to rf. However, as r] increases to of order the Planck mass, the deficit angle 
approaches 2n and static solutions cease to exist ||. Recently, the criterion for topological 
inflation has been made more precise by numerical simulations of the spacetime structure 
around various defects; gauge monopoles ||, domain walls ||, and global monopoles ||[7|]. 
These analyses all focus on defects with unit topological charge, and find that in that case, 
the criterion for topological inflation is 



The aim of this paper is first to verify that the above criterion holds for the important case 
of gauge cosmic strings, and second, to investigate how the criteria for topological inflation 
depend on the topological charge of the defects considered. 

Our motivations are twofold. First, for defects of unit winding, (|1 . 1|) implies that, to 
realize topological inflation, we must work very close to the Planck scale, at which our field 
theories may not be valid. Further, it is clear that topological inflation cannot occur at 
the GUT scale for unit winding defects. We hope that, when higher topological charges 
are included, these constraints will be alleviated and the value of r/ cr will decrease. In fact, 
one might expect such behavior from considering the static string solutions, since the deficit 
angle we mentioned above is also dependent on the winding n, and so static solutions should 
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cease to exist for lower values of r/ if n > 1. The spacetime structure for such higher winding 
strings is the focus of this paper. 

Our second motivation comes from models in which particles are described as solitons 
||, and in particular from the dual standard model |§. In this theory, all the standard model 
particles arise as monopoles of some bosonic field theory. In such a model it is natural to 
wonder what happens to matter at high densities when the core structure of the solitons 
becomes important. For example, stars can be seen as collections of huge numbers (~ 10 57 ) 
of monopoles. If the monopoles become squeezed together tightly enough for a large region 
of the star to be in the false vacuum with high winding, might topological inflation occur? 

In the present work, fueled by the above considerations, we consider the simple example 
of an Abelian-Higgs vortex with winding n in (2 + 1) spacetime dimensions. We do so 
because it is easier to work with multiple winding vortices than the analogous monopoles, 
although the problem of monopoles for n > 1 is under consideration. In the next section we 
present the model and the equations of motion we solve. We give our initial conditions and 
describe how we expect solutions to behave in some asymptotic regimes. In section III, we 
briefly discuss the implementation of the numerical algorithms we use to solve the equations 
and in section IV we present our results. Section V contains a discussion of the results and 
their implications for topological inflation. 

II. THE MODEL 

Consider the Abelian Higgs model with a complex scalar field $ and a U(l) gauge field 
v4 M , coupled to gravity in (2 + 1) spacetime dimensions. The action is 

where 1Z is the Ricci scalar and the Lagrangian density for the matter fields is 

C = (D^)*D^ - ^F^F^ - V($) . (2.2) 

Here, the covariant derivative is = (V M + ieA^Q, the gauge field strength is F^ = 
d fl A u — d u A^, and the scalar potential is 

7(*) = ^(M-ir l ) , (2.3) 

with A, e constants and rj a mass scale. 

In cylindrical polar coordinates (r, 6), we make the usual Nielsen-Olesen string ansatz 
for the fields 

$(x) = f(r,t)e me , 

A e = -[h(r,t)-n] , 
e 

where the integer n is the winding number of the string. The metric ansatz is 

ds 2 = d f _ e B(r,t) dr 2 _ e C(r,t) r 2 d9 2 ? (2 4) 
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where B and C are functions of the radial coordinate and time. 

This ansatze leads to four Einstein equations and two field equations for a total of 
four unknown functions: B, C, f and h. Two of the six equations are first order in time 
derivatives and are the constraint equations. The four equations we solve are: 
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f + ^(f 2 -V 2 )f = 0, (2.7) 



h - e- B h" - ~(C - B)h + — (C + B')ti + e- B - + 2e 2 f 2 h = 
2 2 r 



(2.8) 



where a dot (prime) denotes a derivative with respect to time (r). At all times, we insist 
that the two additional constraint equations 



BC + e- B (-2C" - C 2 + B'C - 4— + 2— 
\ r r 



32tt 



„— C P —B—C 

B f a , e ^ 2+ e 



2e 2 r 2 



2e 2 r 2 



^ /2 + ^ 2 / 2 + t(/ 2 -^ 2 ) 2 



(2.9) 



and 



2 (7'_(7C / + J BC / -2- + 2- = ^ 



771: 



P L 



-C 



Vf + -^hti 



(2.10) 



are satisfied. 

Now consider the initial conditions for these equations. We begin with a cylindrically 
symmetric string configuration for the fields / and h, which is initially static. We define 
the metric to be initially flat but with non-vanishing first time derivatives. For the metric, 
these conditions are simply implemented as 



B(r,0) = C(r,0) = , 



(2.11) 



with B(r, 0) and C(r, 0) obtained from the constraint equation ( |2.9| ). For the fields, however, 
we need the initial profile functions f(r, 0) and h(r, 0). To simplify this process, we work in 
the Bogomolnyi limit defined by 



(2.12) 
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In this case, there are no forces between static strings and the energy saturates a topological 
bound. Since we are primarily concerned with gravitational effects, we do not expect results 
for f3 7^ 1 to be significantly different. In the Bogomolnyi limit, the static field equations 
reduce to the two first order equations 

f + e (B-ofthl = Q (213) 
r 



h' + ^ B+c ^ 2 (f 2 -r ] 2 ) = 0, 



(2.14) 



which can be solved numerically. 

Our procedure is as follows. We first complete our initial conditions by solving (j2.13||TT4D 
subject to the boundary conditions 



lim f(r) =ii 

Mo) = - 



n 



(2.15) 



We then solve equations ( |2.5| . |2.6| . [2~7|J2.8| ) with the initial conditions we have just described. 
Throughout the evolution we verify that the constraint equations ( |2.9| , |2.10D are satisfied at 
each step as a check of our numerical scheme. 

For a given topological charge n, this procedure is performed over a range of values of the 
symmetry breaking scale i]. We define a solution exhibiting topological inflation to be one for 
which the total physical volume, V*(t), in the core of the defect is increasing exponentially. 
We define V*(t) by 



V* (t) = 2n / drr exp 
Jo 

where the core radius, r*(i) is defined by 



B(r,t) + C(r,t) 



(2-16) 



f[n(t)} = 1 . 



(2-17) 



Determining the functional form of r) cr (n) is the central result of this paper. 

It is a useful check of our results that one may simply derive an upper bound for the 
expansion rate in the core. Assume that inflation occurs and that the metric components B 
and C become very large compared to other fields in the core of the defect. Further, assume 
that only the vacuum energy of the scalar field is important in the core, i.e. that = 
with no derivatives important there. In this approximation, the equations of motion for the 
metric simplify dramatically and are easily solved to give 



B(0, t) ~ C(0,t) 




(2-18) 



Thus, any inflationary behavior we observe should have an associated Hubble constant H 
that satisfies 
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H < (8ttA) 1/2 ( (2.19) 
\m p J 

Our intuition for believing that higher topological charges will alleviate the high symme- 
try breaking scales required for topological inflation comes from two sources. First, consider 
the asymptotic form of the metric for static cosmic string solutions (in 2+1) dimensions 

ds 2 2+1 = dt 2 - dr 2 - r 2 d9 2 , (2.20) 

where 9 is the angle in a locally Minkowski but globally conical spatial section, taking values 
in the range 

< 9 < 2tt ( 1 - 4|d-0 . (2.21) 
V m pJ 

For strings of unit winding, this metric is applicable as long as the deficit angle is less than 
27r, that is, for symmetry breaking scales rj <C m p . However, for higher winding, static 
solutions cease to exist for 




(2.22) 



Thus, we expect that asymptotically static solutions become impossible at a lower critical 
symmetry breaking scale for defects with multiple winding and we might guess that the 
critical value of r\ at which static solutions cease to exist falls off as 1/ yfn. It is natural to 
wonder if the same is true in the core of these defects, although, of course, the absence of 
static solutions does not guarantee that the core will inflate. 

Second, a perturbative analysis of the matter fields around the center of such defects 
demonstrates that defects with |n| > 1 have a wider core and higher energy density than 
the corresponding unit charge configurations. Both these effects suggest that topological 
inflation might be more easily achieved in high winding defects. Unfortunately, it does not 
seem possible to quantitatively understand the effect of multiple winding on topological 
inflation with an analytic approach. Thus, here we have solved the system numerically, in 
the spirit of other authors in the unit winding case ]7|,|5]J| . 

III. NUMERICAL IMPLEMENTATION AND RESULTS 

In this section, we present the results of our numerical simulations of the (2 + 1) dimen- 
sional Einstein- Abelian-Higgs system. The system of non-linear partial differential equations 
we study are non-trivial to solve numerically. Therefore, before we present our results, let 
us briefly discuss the numerical techniques we use, in the hope that this discussion will help 
others investigating similar problems. 

A. Numerical Implementation 

There are two stages to solving the equations. The first is to generate the field configura- 
tions that will serve as the initial values for the time dependent evolution equations, and the 
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second is to integrate the partial differential equations that describe the time dependence of 
the fields. To attack the former, we observe that the Bogomolnyi equations ( p. 13 ) and Q2.1-1 ) 
have asymptotic solution 

lim/(r)~/or", (3.1) 
r— »0 

where /o is a constant of integration. This constant can be used as a free parameter in a 
shooting method solution to the boundary value problem. We find that shooting is effective 
in generating accurate solutions out to a radius of lOr], sufficient for both fields to reach 
their r — > 00 asymptotic values for the windings we consider. 

With the initial conditions in hand, we must now consider how to integrate the time 
dependent partial differential equations. Typically one replaces derivatives with finite dif- 
ference approximations. For a generic variable y(r,t), one solves for values on a lattice 
y(r,t) — ► yf, where subscripts indicate the position in the space lattice and superscripts 
indicate the location in the time lattice. Derivatives are replaced with finite difference 
approximations, e.g. 



yf +1 - 2yj + yj 1 



dt 2 



y" 



yl+i - %i + yl-i 



dr 2 



which are second order approximations in dt and dr. However, there is in general no guaran- 
tee that a particular differencing scheme is stable. That is to say, for poor schemes the result 
from integrating the difference equations may diverge exponentially from the true solution. 
To test the differencing methods, we may use the stability analysis for linear equations which 
is covered in any good reference on partial differential equations (see e.g. |TU[ chapter 19, 



and references therein). The results of such an analysis also provide good intuition when 
dealing with the non-linear equations for the metric and fields we are considering. For our 
equations, we find that when solving for / and h, stability is assured if the spatial derivatives 
of these fields are evaluated implicitly. That is, for the j th time step we evaluate the second 
spatial derivative of / as 

f" = fei 1 - 2yl +1 + y£l) /dr 2 , (3.2) 



and similarly for the others. At each time step, the implicit scheme gives us a set of 
coupled linear equations for the fields // + and h{ + , solving which reduces to inverting 
a tri-diagonal matrix, a standard problem in linear algebra. We also find that care is 
needed when evaluating the metric equations ( |2.5|) and ( [2.6|) . Both can be written in the 
form 

C + C 2 = F, (3.3) 

where F just represents the terms on the right hand side of the equation. Treating C as 
an independent variable such that C*/ +1 = C*/ _1 + 2dtC{, we found that it is necessary to 



evaluate C implicitly for eq. (|3.3| ) to be stable, i. e 
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(3.4) 



The above quadratic has two solutions, but one may easily obtain the right one by taking 
the limit as dt — > and noting that, for the correct root, C{ +1 — C\~ should vanish. These 
suggestions worked well for us, although they do not represent the only stable differencing 
schemes and they may not generalize to other similar problems. 

B. Results 

As we mentioned in section II, our strategy was to evolve the initial configurations for a 
given value of the winding n, for various values of the scale rj. As an example, in Figure ([I]) 
we show the metric fields B(r, t) and C(r, t) for an n = 1 string with symmetry breaking scale 
7] = 0.2m p > i] cr . For comparison, in Figure (0) we show the same fields for an n = 5 string 
with rj = 0.07 > r) cr . These plots demonstrate that the metric fields grow exponentially in 
the core but that the core size decreases exponentially. It is the competition between these 
two effects that determines whether inflation occurs or not. 

In both cases inflation is occurring in the core of the defect, although this is not clear until 
we apply our criterion that the total volume of physical space be increasing exponentially 
in the core. Note that, as an artifact of our initial conditions, there is an initial period of 
time during which the system relaxes to its final state. 

To illustrate how the criterion is applied, Figure @ shows V /V as a function of time 
t for two values of r], one for which the core inflates, and the other for which it does not, 
for two cases of the topological winding, n = 1 and n = 5. It is the qualitative difference 
between these two classes of scales that allows us to home in on r] cr ~ 0.16 for n = 1 and 
rjcr ~ 0.06 for n = 5. In the n — 1 case, the lower curve appears to turn up at late times. We 
believe this to be due to our rigid definition of the core radius, which does not take account 
of oscillations that appear in the matter fields. However, finite-size effects in the simulations 
limit our ability to test this. 

Finally, in Figure (^) we show the relationship between rj cr and the topological charge n 
on a log- log plot. The points are best fit by a linear relationship 

r] cr ~ an p , (3.5) 

where 

a = 0.16, p=-0.56. (3.6) 

This is in excellent agreement with the naive estimate p = —0.5 obtained by analyzing the 
point at which the static asymptotic metric ceases to exist. 

IV. CONCLUSIONS AND DISCUSSION 

We have analyzed the onset of topological inflation in the cores of cosmic string solutions 
to the Einstein- Abelian-Higgs system in (2 + 1) dimensions. For a soliton in a given sector of 
topological charge n, inflation occurs in the defect core if the symmetry breaking scale t] is 
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FIG. 1. The metric fields B{r) and C(r) for an n = 1 string with symmetry breaking scale 
rj = 0.2m p > rf cr . The functions are plotted at equal time steps with the higher amplitude curves 
occurring at later times. 
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FIG. 2. The metric fields B{r) and C(r) for an n = 5 string with symmetry breaking scale 
rf = 0.07m p > r] cr . The functions are plotted at equal time steps with the higher amplitude curves 
occurring at later times. 
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FIG. 3. V/V as a function of time t for two values of rj, one for which the core inflates (upper 
curve in each figure) , and the other for which it does not (lower curve in each figure) , for two cases 
of the topological winding, n = 1 and n = 5 
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FIG. 4. A log-log plot of rj cr , in units of m p , versus the topological charge n. The "error bars" 
denote the ranges within which we could numerically bracket i] cr . Their differing sizes reflect our 
initial trial and error guesses for the bracketing. 



greater than a critical value f] C r{n). The functional dependence of r] cr on the winding n was 
determined numerically and was found to be monotonically decreasing roughly, though not 
exactly, as 1/x/n. This result supports the intuition about defects with multiple winding 
gained from the asymptotic metric of static solutions and from perturbative analyses of the 
core fields. If our results can be extrapolated to very large n, and if strings with such high 
winding form in phase transitions, then it is possible that topological inflation could occur 
in GUT scale defects. 

The present analysis is especially relevant in theories in which particles are viewed as 
solitons. In these theories one would expect that, at high densities, the solitonic nature 
of particles would become important. Our results then show that, provided the number of 
particles is large enough, the tightly squeezed state of particles can start inflating. This 
may be relevant for the gravitational collapse of stars since the number of particles in a 
star is of order 10 57 . In the context of the dual standard model, all particles correspond to 
magnetic monopoles and so we would expect the present considerations to apply there also. 
However, in this model, baryon number is not a conserved quantity and it is possible that 
the star evaporates before inflation can set in, much like in the scenarios recently considered 
in 



11,12 



A number of related investigations are suggested by our analysis. First, is it possible that 
the collision of a monopole and an antimonopole can result in an inflating region? Some years 
ago, Farhi, Guth and Guven |l3j considered the possibility of creating a universe in particle 
collisions. Is their (negative) conclusion applicable even in soliton collisions? Second, we 
have only considered strings at critical coupling. For different choices of couplings, the strings 
could attract or repel each other. How do our results depend on the coupling constants? 
Does the instability of higher winding strings to decay into those of lower winding, come 
into play and terminate topological inflation after a certain number of e-folds? We hope to 
return to some of these questions in future investigations. 
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